      SUBROUTINE SETUP
      RETURN
      END
C   HS76
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=X(1)**2+.5D0 *X(2)**2+X(3)**2+.5D0*X(4)**2-X(1)*X(3)+X(3)*X(4)
     1   - X(1)-3.D0*X(2)+X(3)-X(4)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)= X(1)+X(1)-X(3)-1.D0
      GRADF(2)= X(2)-3.D0
      GRADF(3)= X(3)+X(3)-X(1)+X(4)+1.D0
      GRADF(4)= X(3)+X(4)-1.D0
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),B(4,3),C(3),GXI,S
      SAVE B,C
      DATA B /-1.D0,-2.D0,-1.D0,-1.D0,-3.D0,-1.D0,-2.D0,1.D0,0.D0,
     1         1.D0, 4.D0, 0.D0/
      DATA C/  5.0D0,  4.0D0,-1.5D0 /
      IF ( GUNIT(1,I+NH) .EQ. -1 )CRES(I+NH)=CRES(I+NH)+1
      IF(I .GT. 3)      GOTO 200
      S=C(I)
      DO      100      J=1,4
      S=S+B(J,I)*X(J)
 100  CONTINUE
      GXI=S
      RETURN
  200  CONTINUE
      GXI=X(I-3)
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADGI(NX),B(4,3)
      SAVE B
      DATA B /-1.D0,-2.D0,-1.D0,-1.D0,-3.D0,-1.D0,-2.D0,1.D0,0.D0,
     1         1.D0, 4.D0, 0.D0/
      IF ( GUNIT(1,I+NH) .EQ. 1 )RETURN
      CGRES(I+NH)=CGRES(I+NH)+1
      IF(I .GT. 3)    GOTO 200
      DO      100      J=1,4
      GRADGI(J)=B(J,I)
  100 CONTINUE
      RETURN
  200 CONTINUE
      DO      300      J=1,4
      GRADGI(J)=0.D0
  300 CONTINUE
      GRADGI(I-3)= 1.D0
      RETURN
      END
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS76'/
      DATA (X(I),I=1,4)/ 4*.5D0/
      DATA N/ 4/ , NH/0/ , NG/ 7/
      DATA DEL0/1.00D0/ ,TAU0/1.D0/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,7)/-1,0,0,-1,0,0,-1,0,0,-1,0,0,
     F           1,1,1,1,2,1,1,3,1,1,4,1/
      END

